Instabilities in numerical loop quantum cosmology 



Jessica Rosen, Jae-Hun Jung and 

Mathematics Department, University of Massachusetts at Dartmouth, North Dartmouth, Massachusetts 0274^ 



O 
O 



> 

^J- 
O 
l> 

o 

\o 
o 

"o 

cr 
i 

S-H 



X 



Gaurav Khanna 

Physics Department, University of Massachusetts at Dartmouth, North Dartmouth, Massachusetts 0274^\ 

(Dated: February 7, 2008) 

In this article we perform von Neumann analysis of the difference equations that arise as a result 
of loop quantum gravity being applied to models of cosmology and black holes. In particular, 
we study the numerical stability of Bianchi I LRS (symmetric and non-symmetric constraint) and 
Schwarzschild interior (symmetric constraint) models, and find that there exist domains over which 
there are instabilities, generically. We also present explicit evolutions of wave-packets in these 
models and clearly demonstrate the presence of these instabilities. 
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I. INTRODUCTION 

Loop quantum cosmology (LQC) is a symmetry- 
reduced application of loop quantum gravity [2|, a the- 
ory that leads to space-time that is discrete at Planck 
scales. An important result of LQC is that it is free of 
singularities. This is seen as a general consequence of the 
quantum evolution equation, which is a "difference equa- 
tion" for the wave function and does not break down 
where the classical singularity would be . In addition, 
LQC has led to the emergence of a possible mechanism 
for inflation yL development of effective semi-classical 
Hamiltonians |5| and several other interesting directions 
for research. 

Solving these difference equations that arise in LQC 
models is an active area of research today. Notions 
of semi-classicality and the related concept of pre- 
classicality play a very significant role currently, because 
in the appropriate limits, solutions to these LQC models 
must yield the known and expected classical solutions. 
Recently there has been a lot of progress on this front 
and several techniques of both analytical @ and numer- 
ical 0@ nature have been developed to obtain appro- 
priate solutions to LQC models. In this article, it is not 
our focus to further develop methods for obtaining solu- 
tions to LQC models. Instead, we investigate a specific 
property, i.e. the "generic stability" of these solutions in 
certain homogeneous LQC models. Note that a study of a 
similar type, but in a more general context and emphasis 
is available in literature @- In this article, we concen- 
trate on Bianchi I locally rotationally symmetric (LRS) 
cosmology and the Schwarzschild interior geometry. 

While difference equations can be solved very natu- 
rally by numerical simulations on computers, it is essen- 
tial to note one crucial difference between LQC models 
and the usual finite-difference approaches towards solving 
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and analyzing common differential equations. An impor- 
tant aspect of the numerical analysis of differential equa- 
tions is the flexibility of choosing various finite-difference 
schemes, that in the continuum limit yield the differ- 
ential equation under consideration. In fact, it is this 
freedom that is used to develop sophisticated algorithms 
for numerical evolutions that are stable and convergent. 
Unfortunately, in numerical LQC we lack this freedom. 
This is because it is the difference equation itself that is 
fundamental, and is not the result of a discretization of 
some differential equation. However, there is some lim- 
ited flexibility, that is related to quantum ambiguities 
and different variations of the Hamiltonian constraint in 
the theory. In this article, we investigate two different 
forms of the Hamiltonian constraint for Bianchi I LRS 
model and the dependence on a quantum ambiguity in 
the Schwarzschild interior geometry model. 

This article is organized as follows. We first present 
the von Neumann stability analysis of the non-symmetric 
constraint in the context of the Bianchi I LRS LQC 
model. Then, we study the same for the Schwarzschild in- 
terior LQC model using the symmetric Hamiltonian con- 
straint. We also make some remarks about the Bianchi I 
LRS symmetric constraint, which is a special case of the 
Schwarzschild model. Finally, we end with a discussion 
of our results. 



II. NUMERICAL STABILITY ANALYSIS OF 
LQC MODELS 

A. Bianchi I LRS non-symmetric constraint 

The Bianchi I LRS homogeneous model has two de- 
grees of freedom which are given by two connection 
components (A, c) and the conjugate momenta (pa,Pc) 
with symplectic structure given by {^4,pyi} = ^jk and 
{c,p c } — jk. Here, n = 8ttG is the gravitational constant 
and 7 the Barbero-Immirzi parameter of loop quantum 
gravity. The momenta (j>a,Pc) are components of a den- 
sitized triad which determine the scale factors ai of a 
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a 2 , a 3 = 
■dy 2 )+p 2 A /\Pc\dz 2 



Bianchi I metric by a\ — yj\p c \ = 
Thus, the metric is ds 2 = \p c \(dx 2 
Cartesian coordinates. 

After quantizin g th e model with techniques from loop 
quantum gravity [lOj, the triad components pa and p c 
become basic operators with discrete spectra, pA\m, n) — 
\^ipm\m,n) and p c \m,n) = ^Ppr^m^n). The wave 
function is thus supported on a discrete minsuperspace, 
and is a solution to the quantized Hamiltonian constraint 
(non-symmetric) subject to a two-parameter difference 
equation, 
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in matrix form as seen below, 
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verse matrix of the LHS matrix yields, 
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d 2 (m) = — , and d 2 (0) = 0, 
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and d(0) = 0. 



We shall interpret the discrete parameter, n as a quan- 
tum analog of physical "time" and the parameter, to as 
"space". We shall focus on the numerical stability of 
the given difference equation We employ von Neu- 
mann analysis which is used to investigate the stability 
of generic finite difference schemes. The procedure is 
to perform a (discrete) Fourier transform along all spa- 
tial dimensions, thus reducing the system to a recursion 
relation in time. Using the relation, one computes the 
"amplification" or "gain" by calculating the ratio of the 
values of the transform at consecutive time-steps. If this 
amplification ratio's absolute value is bounded by unity, 
the system is stable . We begin by writing the above 
equation in the following fashion. Let, 
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then the equation can be expressed in the form, 
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Now, as prescribed for standard von Neumann stability 
analysis, we assume that = v p exp(imcj) and = 
w p exp(zTOw), and then the difference equations for vf n 
and are given by 
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Let the coefficient matrix in the RHS of the equation 
above be M(a>,p, fi). Then (v p ,w p ) T becomes 



= JjM(w,i,m) 

3=1 



W 



The necessary condition for stability is that absolute val- 
ues of the eigenvalues of this matrix are bounded by 
unity. These eigenvalues Ai^ of M are given by 

X——D(p, m) sin 2 oj± — \J sin 2 uiD(p, m) [D(p, to) sin 2 lu — 4] , 
where 



D(n, to) 



d(2p)d(2p - I) 
d\(m) 



Now, consider the case, D(n, m) sin 2 lj — 4 < 0. In this 
case, is it clear that the absolute value of both eigenvalues 
presented above, is exactly unity. However, if we consider 
the complementary case, i.e. D(n, to) sin 2 lu — 4 > then 
there do exist eigenvalues with absolute values that are 
larger than unity, which indicates the presence of an in- 
stability in this discrete system. It is interesting to note 
that it is the larger values of to that induce the insta- 
bility. In particular, the instability arises if to is in the 
interval given by (we have converted back to the index 
parameter, n here) 



to 2 > 



sin iod(n)d(n — 1) 



(2) 



where sin to ^ 0. If sin lu = 0, we have immediately 
A = 1. To see the domain of the instability more clearly, 
in figure I we plot the unstable eigenvalue as a function 
of n and to. Here we simply chose sin lo = I. The domain 
of instability is roughly characterized by to > An as seen 
clearly in the figure and also from the inequality that 
appears above. 
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FIG. 1: Unstable eigenvalue of the "evolution matrix" for the 
Bianchi I LRS LQC model with the non- symmetric Hamilto- 
nian constraint. Notice the range where the eigenvalues are 
greater than the unity. 



B. Schwarzschild interior model with symmetric 
Hamiltonian constraint 

Recently, the interior of Schwarzschild geometry was 
quantized using techniques from LQC This is pos- 
sible because the interior space-time of a Schwarzschild 
black hole can be considered as a homogeneous cosmol- 
ogy model of Kantowski-Sachs type. The symmetric 
Hamiltonian constraint in this case, takes the form (for 
r,n>l) 

Mr) KX\ %±{) Mr) K+l - %Z\) 

+ di(r) KM*;+2-™iM*; + moW*;-2) =0 

where 
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m (/i) = A* - 1 



mx(n) = |U 

and k is a parameter that is related to the Immirzi con- 
stant and a quantum ambiguity In this model, we 
shall interpret the parameter r as "quantum time" and 
fi as "quantum space" . Now, let us proceed with von 
Neumann stability analysis of this difference equation. 
Let 

= * r exp(io;/i). 

Then the constraint becomes 

2isin(w)d 2 (T)* r+1 - 2isin(w)d (r)* r " 1 
+ di (t) (h(2 cos(2cj) -K) + 2i sin(2w))) * r = 

As before, we define the auxiliary variables E and B such 
as 

t = ---,-l, 1,3,-.-, 

and 

B l :=¥, t= •••,0,2,4,-... 

We also define the coefficient variables G2(u>, r), Gq(u>, t) 
and H[u>, fi) such as 

GzioJjT) :— 2isbx(u>)d2(T), 
G (uj,t) := 2isin(w)d (r), 

and 

H{u,n) := di(r)(^(2cos(2w) - k) + 2isin(2o;)). 

Then the difference equation can be rewritten in a matrix 
representation 

/G 2 (o;,t) H(u,n) \(ET+i\ 
^ G 2 (w,t-1) )\ ) 



( Go(«,r) \/£ T -'\ 

^-JT(w,aO G (w,r-1) J y B T_2 y 

Let the coefficient matrix M(u>, r, ^t) be 

(G 2 (u,t) ff^/x) rVGo(w,T) 

^ G 2 (u,t-1) J \-H(w,ii) G (w,r-l) y 

Then we have 

(^ 1 )=n M (^)(£) 

The eigenvalue expressions of this matrix are too com- 
plicated to present in this article, and will not be es- 
pecially useful. However, again much like the Bianchi I 
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LRS non-symmetric constraint we discussed before, we 
find that there are regimes over which the absolute val- 
ues of the eigenvalues exceed unity and thus indicate the 
presence of an instability. This happens generically, for 
any value of k. We present below, in figures 2 and 3, 
plots of the unstable eigenvalue as a function of t and 
[a and for different values of k. Note that the k = 1 
case corresponds to the Bianchi I LRS LQC model with 
the symmetric constraint. The domain of instability is 
again roughly characterized by [i > At as seen clearly in 
the figure. The plots look very similar for other values 
of k, so to save space, we do not include more of those 
but point out that with larger values of k the unstable 
region gets progressively larger and the instability gets 
worse (see figure 4). For the sake of completeness, we 
also show a plot of the stable eigenvalue for the case of 
k = 2 in figure 5. 




FIG. 2: Unstable eigenvalue of the "evolution matrix" for the 
Bianchi I LRS LQC model with the symmetric Hamiltonian 
constraint. Notice the range where the eigenvalues are greater 
than the unity. 



C. Evolution of wave- packets 

Below in figure 7 and 8, we also present explicit nu- 
merical evolutions of wave-packets for the k = 2 case. 
These evolutions were performed on the basis of the sten- 
cil depicted in figure 6, which shows the values of the 
wave function that are related by the Hamiltonian con- 
straint in consideration here. Also note that there are 
some residual gauge symmetries that indicate that the 
wave function is unaffected by the gauge transformation 
[i — » — (i, coming from the Gauss constraint Putting 
this together, we can fix the boundaries of our numerical 
simulation to be u = ±M, for large M. Further, there 
are arguments [Tj for choosing W — > as /i — > oo. Thus, 
as long as the maximum grid size M is large enough to 




FIG. 3: Unstable eigenvalue of the "evolution matrix" for the 
Schwarzschild interior LQC model with the symmetric Hamil- 
tonian constraint (k — 2). Notice the range where the eigen- 
values are greater than the unity. 
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FIG. 4: Eigenvalue plotted as a function of k at fixed values 
of /j, — t — 100. Clearly, larger values of k make the sys- 
tem more unstable and also make the unstable region become 
progressively larger. 



avoid sizeable errors in solving the difference equation, 
we can set ^^. M = 0. 

Our results are as expected. If the wave-packet is 
evolved over a domain where the evolution is expected 
to be unstable, we clearly see that the numerical evolu- 
tion fails after a few time-steps. On the other hand, over 
the stable domain we appear to be able to numerically 
evolve the wave-packets for a long time with success. 



FIG. 5: Stable eigenvalue of the "evolution matrix" for the 
Schwarzschild interior LQC model with the symmetric Hamil- 
tonian constraint (k — 2). 



FIG. 6: Stencil, that shows the values of the wave function 
that are related by the Schwarzschild interior or Bianchi I 
LRS symmetric Hamiltonian constraint. 



III. DISCUSSION 



FIG. 7: Unstable wave-packet numerical evolution for the 
Schwarzschild interior LQC model with the symmetric Hamil- 
tonian constraint (k = 2). We chose the range of /i to 
( — 1000, 1000) and we evolved forward in r from 150 to 186. 
The initial wave-packets were centered at ±600 and had a 
width of 50 which places part of them in the unstable region. 




In this article, we performed von Neumann numeri- 
cal stability analysis of some of the difference equations 
that arise in homogeneous LQC models. In particu- 
lar, we focussed on Bianchi I LRS cosmology and the 
Schwarzschild interior geometry. We investigated two dif- 
ferent forms of the Hamiltonian constraint for Bianchi I 
LRS model and the dependence on a quantum ambiguity 
in the Schwarzschild interior geometry model. In all these 
cases, we discovered that there arc large regions in space- 
time that have instabilities, generically. These instabili- 
ties are characterized by rapid exponential growth, and 
therefore they are difficult to interpret physically. The 
presence of these instabilities, especially in the large vol- 
ume limit is troubling because it is indicative of the fact 
that the models do not have the correct semi-classical 
behavior. 

These results may be related to another recent trou- 



FIG. 8: Stable wave-packet evolution for the Schwarzschild in- 
terior LQC model with the symmetric Hamiltonian constraint 
(k — 2). We chose the range of fi to ( — 1000, 1000) and we 
evolved forward in r from 200 to 300. The initial wave-packets 
were centered at ±600 and had a width of 50 which places them 
in the stable region. 



bling result from the study of pre-classical solutions of 
the Schwarzschild interior LQC model using generating 
function techniques |l3j . There, an upper bound was 
discovered on the quantity wc called k here, which in 
turn led to a bound on the value of the Immirzi param- 
eter. This bound on the Immirzi parameter appears to 
be violated by its currently accepted values. It is possi- 
ble that our current results are simply further indicating 
that these models appear to be lacking in some fashion. 
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Another possibility is that these regions that appear 
to have the instability are "forbidden" regions in some 
manner. What is interesting about this possibility is that 
these regions are forbidden quantum mechanically, but 
not forbidden classically (see classical trajectories of this 
LQC model in 11]). If this would be true, it would be a 
very unusual scenario. We do not feel that this is a very 
likely possibility. 

Finally, we'd like to comment on yet another interest- 
ing possibility 0]. Note that a region is semi-classical 
if the length and time scales are large and the curva- 
ture scales are small. In mini-superspace models, such as 
the ones under consideration here, the relevant curvature 
can become large with the spatial size. This creates the 
kind of instability problems we're observing here. More- 
over, difference equations that are derived in the context 
of mini-superspace models do not take into account the 
subdivision of discrete geometry which is expected when 



the Universe grows larger. Taking this into account re- 
quires introducing inhomogeneity or introducing a scale- 
dependent discretization, which is currently a topic of 
exploration. However, even in these cases its is possi- 
ble that unstable regions may not go away completely 
and thus may still present interesting opportunities for 
further research. 
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